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Abstract 

The NOAA active region AR 11029 was a small but highly active sunspot region virhich 
produced 73 GOESsoit X-ray flares during its transit of the disk in late October 2009. 
The flares appear to show a departure from the well known power-law frequency-size 
distribution. Specifically, too few GOES C-class and no M-class flares were observed 
by comparison with a power-law distribution (Wheatland in Astrophys. J. 710. 1324, 
2010). This was conjectured to be due to the region having insufficient magnetic energy 
to power the missing large events. We construct nonlinear force-free extrapolations of 
the coronal magnetic field of active region AR 11029 using data taken on 24 Octo- 
ber by the SOLIS Vector-SpectroMagnetograph (SOLIS/VSM), and data taken on 
27 October by the Hinode Solar Optical Telescope SpectroPolarimeter (Hinode/SP). 
Force- free modeling with photospheric magnetogram data encounters problems because 
the magnetogram data are inconsistent with a force-free model. We employ a recently 
developed 'self-consistency' procedure which addresses this problem and accommodates 
uncertainties in the boundary data (Wheatland and Regnier in Astrophys. J. 700, L88, 
2009). We calculate the total energy and free energy of the self-consistent solution which 
provides a model for the coronal magnetic field of the active region. The free energy 
of the region is found to be w 4 x 10^^ erg on 24 October, and « 7 x lO'^^ erg on 
27 October. An order of magnitude scaling between RHESSI non-thermal energy and 
GOES peak X-ray flux is established from a sample of flares from the literature and 
is used to estimate flare energies from observed GOES peak X-ray flux. Based on the 
scaling, we conclude that the estimated free energy of AR 11029 on 27 October when 
the flaring rate peaked is sufficient to power M-class or X-class flares, and hence the 
modeling does not appear to support the hypothesis that the absence of large flares is 
due to the region having limited energy. 

Keywords: Active Regions, Magnetic Fields, Corona 



1. Introduction 



The dynamics of the Sun's diffuse outer atmosphere, the corona, is dominated by the 
solar magnetic field (jAschwandenl . 120041 ). Solar flares are the result of the explosive re- 
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Figure 1. Evolution of active region AR 11029 from 23 to 27 October 2009. The figure shows a 
collage of magnetograms from the Michelson Doppler Interferometer (MDI) aboard the Solar and 
Heliospheric Observatory. The data show the linc-of-sight magnetic field 120 km above the base of 
the photosphere, and the field values saturate at ±500 gauss (G). The times when the data arc taken 
are indicated in the panels. 



lease of magnetic energy from intense coronal magnetic fields associated with sunspots. 
Understanding the storage and release of magnetic energy in the corona is fundamental 
to understanding solar flares, and other coronal transient events {e.g. coronal mass 
ejections and transient X-ray brightenings). 

The NOAA solar active region AR 11029 emerged on 21-22 October 2009, began 
flaring on 24 October 2009, and rotated off the disk on 1 November 2009. Figure [T] 
shows the evolution of the photospheric magnetic field of AR 11029 between 23 and 
27 October as observed by the Michelson Doppler Interferometer (MDI) aboard the 
Solar and Heliospheric Observatory. Over the period shown the region grew in size and 
complexity. A detailed summ ary of the region's evolution during its transit of the disk 
is given bv lWheatlandl ^01^. 



Active region AR 11029 was highly active, pr oducing a large number of sm all flares 
during a period of pronounced solar minimum ( Livingston and Penn . 20091 ). In total 



73 soft X-ray events are listed in the event catalog]] compiled from observations by the 
Geostationary Operational Environmental Satellites (GOES) during the region's transit 
of the disk. The majority are small GOES B-class events (peak 1-8 A soft X-ray flux in 
the range lO'^Wm-^ to lO'^Wm'^). The largest is a GOES C2.2 event (peak soft X- 
ray flux of 2.2 x 10~^Wm~^). No moderate sized M-class flares (peak flux 10^^ Wm~^ 
to 10""' Wm^^) or large X-class flares (peak flux > lO^'^Wm"^) were produced. 

The frequency and size of solar flares (by 'size' we mean a measure of the magnitude 
e.Q. GOES peak flu x) are related by a well known power-law distribution {e.g. Akabanj 



19561 : lHudsoi]ll99lh . This distribution is universal in the sense that the same power law 



index is observed at all times, for events occurring over the entire Sun and for events 



^ The GOES event catalog is found online at |http://www.ngdc.noaa.gov| 
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that occur in different active regions (|Wheatlandl . l2000[) . The power law in s ize reflects 
an underlying power law in the distribution of flare energies (jHudsonl . Il99lh . 

Analysis of the statistics of the flares produced by active region AR 11029 using 
the GOES event data provides strong evidence for a depa rture from the exp ected 
power-law frequency-peak flux distribution for this region (IWheatlandl . 120101). The 



region produced fewer large events than expected given the number of small events, 
based on the power law distribution, and hence did not appear to follow the universal 
distribution. 

The upper panel of Figure [5] shows the frequency-peak flux distribution for the flares 
produced by AR 11029 as a cumulative histogram (diamonds), and the curves show two 
model distributions. The straight line corresponds to a power-law model distribution 
P{S) ^ where S is the peak flux and 7 is the power-law index. The curved line 
corresponds to a power law model with an exponential rollover at large event sizes 
P{S) ^ S'^'*' cxp(— S'/cr) where A is the power-law index for small S and a is the 
rollover size. The vertical dotted line shows the peak flux above which the two models 
are assumed to apply (10~''Wm~^). There is a noticeable departure from the power 
law at the large end of the scale. Given the data, a Bayesian model comparison taking 
into account all events showed that the rollov er model was favo red over the simple 
power-law model with an odds raticH of « 220 jWheatlandl. IgOlOl) . 

The absence of large flares was attributed bv I Wh eatland (j2010[) to the region having 
insufficient energy to produce large events. In principle an upper bound on flare size 
must exi s t whi ch reflects the finite energy stored in the active region's magnetic field 
(jHudsonl . Il99l[ ). However, it has proven difficult to observationally identify this effect 
in flare statistics (jHudsonl . l2007l ) . 

By determining the energy stored in the coronal magnetic field of active region AR 
11029 it is in principle possible to test the hypothesis that the absence of large flares 
results from a lack of energy stored in the active region's magnetic field. The energy 
available for release in flares is called the free energy E't, and is stored in the component 
of the coronal magnetic field generated by large scale electric currents in the corona 
(jPriest and Forbed . [200I) . Equivalently, the free energy is the energy in excess of the 
current free (potential) component i.e. 



El = E - Eo ^ ^ f \B\^dV - ^ I IBopdF, 
2/^0 J 2^0 



(1) 



where E is the total energy of the active region magnetic field B, Eq is the energy 
of the current free (potential) field Bq, and the integrals are over the coronal volume 
of the active region. The potential field is uniquely defined in a finite volume by the 
normal component of the magnetic field on the closed surface surrounding this volume. 

Determining the energy of an active region from Equation (jl} requires knowledge of 
the magnetic field over the three dimensional coronal volume of the region. This infor- 
mation is difflcult to obtain observationally. Solar magne tic fields affect the polarization 
of the Sun's optical emissions through the Zeeman effect ( Landi Deel'Innocenti and Landolfil .[ 
20041 ). This provides one diagnostic tool for inferring the solar magnetic field from 
spectropolarimetric measurements of magnetically sensitive lines, if such lines are avail- 
able. Vector magnetograms are two-dimensional maps of the vector magnetic field 



^The odds ratio is the relative probability that the data favors the rollover model over the power law 
model: a value of one indicates both models are equally li kely, and a va lue of 220 indicates that the 
rollover model is ~ 220 times more likely given the data l lJavned [2003) . 
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Figure 2. The u pper panel shows the frequency-size distribution of soft X-ray flares for active region 
AR 11029 jWhea tland. 201i|), plotted as a cumulative number distribution for the background sub- 
tracted peak soft X-ray flux of G0-B5 events. The straight line shows a simple power law distribution, 
and the curved line shows a power law with an exponential rollover at large peak fluxes. The dotted 
vertical line shows the size threshold above which the models are assumed to apply (10~^Wm~^). 
A Bayesi an model comparis on favors the rollover model to the straight power law with an odds ratio 
of fa 220 jWheatlandl . l20Toh . The lower panel shows the cumulative number of GOES flares observed 
in active region AR 11029 versus time, between 24 October when the region began flaring, and 1 
November when it rotated off the disk. The times when SOLIS/VSM and Hinode/SP magnetogram 
data are available are indicated in the lower panel by the solid vertical lines. 



B close to the level of the photosphere, d erived from spectropolarimctric measure- 
ments (jLandi Degl'Innocenti and Landolfil . 120041 ) . The field values are not coronal, 
but photospheric values. The coron al m agnetic field cannot be determined observa- 
tionally except in special cases (e.g. iHvd cr 196-^, although a nu mber of new meth- 



ods based on radio (e.g. IWhite and Kundijll997t iLee et aZ.lll998l ) and infrared {e.g 



Lin. Kuhn. and Cotilterj[200l observations are under development 



The lack of observational methods for determining the coronal magnetic field moti- 
vates mod eling of the cor onal field from th e photospheric vecto r magnetogram data [see 
reviews bv lSakurail (|l989f ) or more recentlv lwiegelmannl (|2008l ). for an overview of the so 
called 'reconstruction problem'] , and such models permit estimates of the free energy of 
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the c oronal magnetic field fe.q. lRegnier and Pries tl2007l: iThalmann. Wiegelmann. and RaouafiB 



2008f ). The magnetic field on the lower boundary of the coronal volume is provided by 



the spectropolarimetric observations and this formally defi nes a boundary value prob- 
lem f or determining the coronal field, given a suitable model ( Amari. Boulmezaoud. Alvl .[ 



20061 ) ■ In this way the three dimensional coronal magnetic field can be 'extrapolated' 



from two dimensional magnetogram data. 

Nonlinear force-free (NLF F) models of t he coronal magnetic field assume the field 
is static and force free {e.g. ISturroclj 1994 ). This means all non- magnetic forces are 
neglected and the magnetic Lorentz force is assumed to be self balancing, that is 



J X B 0, 



(2) 



where J is the coronal electric current density. In practice a numerical approach is used 
to solve the force-fr ee equations. A number of nume rical methods have be en developed 
(Wiegelmann, 2008 ) and tested on known solutions ( Schriiver et al . 20061) . Section [^T] 
provides further details of force-free modeling. 

For ce- free modelin g enco unters a number of problems in application to magnetogram 
data ( De Rosa et al . 2009f ). This appears to be primarily a result of the inconsis- 
tency between the forc e-free model s and the boundary data. Whilst the corona may 
be close to forc e-free (IGarvl. 120011) , the photosphere is subje c t to significant non- 
magnetic forces (jMetcalf et adll995t[Garvll200l[ ). [Metcalf et a/.l (| 19951 ) estimated that 
the solar atmosphere d oes not become force fr e e unt il at least 400 km above the 
base of the photosphere. IWheatland and Regnier liooi) developed a procedure based 
on an implementation o f the Grad-Rubin method for calculating nonlinear force-free 
fields l| Wheatland . 120071 ) which addresses the inconsistency between the data and the 
model. This approach identifies a force-free magnetic field solution with boundary 
conditions close to, but not exactly matching, the magnetogram boundary condi- 
tions. The resulting 'self-consistent' field may represent the coronal magnetic field. 
The self-consistency procedure can accommodate uncertainties in the boundary data 
(details of this method are given in Section [2.2(1 . This method has previously been 
appli ed to vector magnetogram dat a without incorporating the uncertainties in the 



data (IWheatland and Regn: 



( Wheatland and Lekal . 120111) 



In this paper we use the 



er 



20091 ). and has recently been applied with uncertainties 



Wheatland and Regnieil (|2009l) self-consistency procedure 



to reconstruct the coronal magnetic field of AR 11029 from vector magnetogram data. 
The magnetograms are derived from spectropolarimetric data provided by the US 
National Solar Observatory's (NSO) Synoptic Optical Long term Investigations of 
the Sun Vector-SpcctroMagnetograph (SOLIS/VSM), and the Solar Optical Telescope 
SpectroPolarimetcr onboard the Hinode spacecraft (Hinode/SP). Based on the recon- 
structed coronal magnetic field we calculate the total energy stored in the active region 
and the free energy Ef, defined by Equation ([1]), and compare these results to a lower 
bound on the energy corresponding to the size of the largest observed GOES event 
produced by AR 11029. 

This paper is structured as follows. Section 2 provides background on nonlinear 
force-free models and the self-consistency method, and Section 3 explains the use of 
the SOLIS/VSM and Hinode/SP data. Section 4 explains the application of the force- 
free model to the SOLIS/VSM and Hinode/SP data and presents the results of the 
modeling. Section 5 discusses the results, and Section 6 presents the conclusions. The 
Appendix contains a glossary of term used throughout the paper. 
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2. Background 

2.1. Force- free Modeling of the Coronal Magnetic Field 

The difficulties associated with observational determination of the coronal magnetic 
field have motivated numerical extrapolation methods which aim to 'reconstruct' the 
coronal magnetic field from photospheric spectropolarimetric data. This requires a 
model of the coronal magnetic field. 

The nonlinear force-free model is a static model which assumes non-magnetic forces 
in the corona {e.g. gravity and pressure forces) are negligible compared to the magnetic 
(Lorentz) force. The latter assumption is approximately met when the plasma beta, 
defined as /3 = 2yLtop/|Bp, where p is the gas pressure, is much less than unity. A 
force-free field may be defined as the solutioro to the equations 



B • Va 



(3) 



and 



V X B = aB 



(4) 



(jSturrockl 119941 ). Here B is the magnetic field vector and a is the force- free parameter, 
which is related to the electric current density J = /x^^V x B by the relationship 
J = aB/zip. 

Equations ^ and ^ are nonlinear in the general case (when a is a function of 
position, i.e. a = a(r) where r is the position vector). The nonlinearity, and the mixed 
elliptic-hyperbolic character of the equations makes them difficult to solve, motivating 
a numerical approach. Several numerical methods fo r solving Equat ions ( Bt and dH) in 
a sola r context have been developed [see reviews by Sakurai ( 1989[ ) and Wiegelmamj 
(|2008l )]. A version of the current-fi eld iteration o r Grad-Rubin app roach, originally 
proposed by Grad and Rubin ( 1958 ). is used here ( Wheatlandl . 2007 ). 

Th e bou ndary conditions f or solv ing the force- free equations are described bv lBineaul 
(jl972l ) and lGrad and RubirJ (119581 ). In the solar context, the boundary conditions are 
Bn, the component of the magnetic field normal to the photosphere and ag, the value 
of a on the photosphere over one polarity of the field i.e. over the region where i?„ > 0, 
or the region where Bn < 0. We refer to the two choices, defined by the two different 
polarities, as the P and N choices respectively. 

On the scale of an acti ve region the curvature of the photosphere is often ignored 
( Garv and Hagvard . 1990[ ). and Equations ([3]) and @ are solved in the Cartesian half 
space z > 0, where z = is the photosphere. In this case B„ = B^, and the force- free 
parameter a on the photosphere is defined by 



«o = 



1 



fdBy dB^ 



\ dx 



dy 



(5) 



z=0 



The right hand side of Equation ([5]) may be estimated from photospheric vector 
niagnctogram data, which we assume provides the vector components of B at z = 0. 



^In this paper when we refer to a 'model' we mean the nonUnear forcc-frec model defined by these 
equations. By 'solution' we mean a solution to these equations, for specific boundary conditions. The 
Appendix provides a glossary of terms. 



SOLA: main.tex; 21 October 2011; 0:27; p. 6 



7 



The boundary value problem outlined in this section is for an infinite half space. 
To compute a numerical solution it is necessary to solve the equations on a finite 
grid of points. This introduces the problem that boundary conditions on B and a are 
unknown on the top and sides of the finite domain and force-free extrapolation codes 
must deal with this problem. The code used here assumes the magnetic field is periodic 
in the x and y direction s, and that the force-free parameter a is zero on the top and 
side boundaries ( Wheatland. 2007i) . Equation © requires that a is constant along field 
lines, and hence field lines which leave the top or side of the box have a set to zero along 
their length. This approach deals with the missing boundary data on a and ensures 
that V • J = is satisfied in the volume. The handling of missing boundary data is a 
source of uncertainty for the force-free modeling performed here, and in particular we 
expect the free energy to be decreased by the removal of currents implied by setting a 
to zero. In addition, images of the solution introduced by the periodicity will affect the 
solution. The effect of this should be small when the computational domain is large so 
the region of interest is isolated from the boundaries. 



2.2. A Self-consistent Model of the Coronal Magnetic Field 



Force - free modeling encoun ters problems when applied to photospheric magnetogram 
data ( De Rosa et ad [20091 ). Firstly, the iterative methods used to solve Equations ^ 
and (HI) fail to strictly conv erge if large electric currents are present in the boundary 
data ( De Rosa et al I l2009l) . Secondly, methods using the formally correct boundary 
conditions given in Section 12.11 produce P and N solutions which differ both in terms 
of energy and field line connectivity. The second problem makes it difficult to estimate 
the energy of the region from the model because two solutions with different energy 
are obtained from the same magnetogram. These problems are in part due to noise in 
the data and uncertainties associated with the inversion of polarization measurements 
to give boundary magnetic field values, errors in the resolution of the 180-degrec am- 
biguity in the direction of the field transverse to the line-of-sight (see Section [3TT|) . and 
errors introduced by the process of estimating photospheric values of a. However the 
fundamental proble m is likely the incon sistency between the force-free model and the 



magnetogram data (jMetcalf et all 119951 ) 



The self-consistency procedure developed by IWheatland and Regnier ( 20091 ) ad- 
dresses the inconsistency between the boundary data and the force-free model. The 
method iteratively modifies the values of a on the boundary, which we denote ao, 
until the P and N solutions agree. The self-consistency procedure can also take into 
account the uncertainties on uq which we denote cr. We refer to repeat applications of 

the method involving modifications of ao as self-consistency 'cycles\ 

Here we give a brief outline of the s elf-consistency procedure (see lWheatland and Reeniei^ 



( 2009f ) or Wheatland and Leka (2011) for a detailed explanation). A single self-consistenc};^ 



cycle can be broken into three steps. The first step involves constructing the P and N 
solutions from the the boundary data Bz{x, y) and Qo (x,?/). The Grad-Rubin method, 
an iterative method for solving Equations ^ a nd dH) (iGrad and Rubinl . 1958 ) is used, 
with an implementation in code described by Iwheatlandl ( 2007 ). In general the P 
and N solutions obtained differ indicating the boundary data are inconsistent with a 
force-free model. 

The second step of a self-consistency cycle involves constructing a second set of 
a values on the boundary, which we denote ai, based on the connectivity of the P 
and N solutions. For a force-free field the parameter a is constant along field lines [an 
implication of Equation ([3])]. Closed magnetic field lines connect points in the boundary 
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with i?^ > to points with Bz < 0, which defines a mapping between regions of P 
polarity and N polarity. The mapping defined by the field lines of the P solution maps 
values of ag in the P polarity to the N polarity (since a is constant along field lines). 
This defines a second set of a values over the N polarity. Similarly the N solution 
maps values of ao from the N polarity to the P polarity, which defines a second set 
of a values over the P polarity. Combined the values of ao mapped by the P and N 
solutions cover the entire lower boundary and we denote this second set of values ai. 
In general ao(a;, y) ^ ai{x, y) unless the boundary data arc consistent with a force-free 
model. Uncertainties on ai, which we denote ci, are constructed from a in the same 
way that ai is constructed from ao- 

The third step involves constructing a new set of a values on the lower boundary 
based on ag and ai. The self-consistency procedure treats ao and ai as two separate 
observations of a on the lower boundary and Bayesian probability is used to decide 
on the most prob able value at each boundary point given the uncertainties a and 
CTi ( Javnesl 2003 ). Assuming the uncertainties a (and by construction also cti) are 
Gaussian uncertainties, the new distribution of a on the lower boundary is 



1/0-2 + 1/al 



(6) 



If all uncertainties are equal then anew is the average of the two values at each point i.e. 
ttnow ~ i<^o + cn)/2. An important property of Equation ^ is that anew is unchanged 
under a global rescaling of the uncertainties a and cti, i.e. Equation ([6]) is unchanged 
by the replacement a Ca, where C is a constant. Hence the absolute value of a at 
a point is not important: only the value with respect to other points is important. 

The process of constructing P and N solutions and then applying Equation © 
(all three steps) is called a 'self-consistency cycle'. This procedure is iterated until 
ao = ai at all points ii i the boundary: a single force-f ree solution (a 'self-consistent' 
solution) is obtained. In Wheatland and Regnierl ( 2009t ) the self-consistency procedure 
was demonstrated in application to data from AR 10953, with u ncertainties assumed 
equal at all boundary points, and in IWheatland and Lekal (|201l[ ) the calculation was 
repeated with uncertai nties derived from the s pectro polarimetric inversion. This paper 
follows the approach of I Wheatland and Lekal (j201lr ). 



3. Data and Modeling 

Vector magnetogram data for active region AR 11029 arc available from the ground 
based SOLIS Vector-SpcctroMagnetograph (SOLIS/VSM), and from the Solar Optical 
Telescope SpectroPolarimeter on board the Hinode satellite (Hinode/SP). One SO- 
LIS/VSM magnetogram per day is available for 24, 25, 26 and 28 October, and one 
Hinode/SP magnetogram per day is available from 24-29 October. All magnetograms 
are derived from observations of the Fe i (630nm) multiplet which provides the vector 
magnetic field close to the level of the photosphere (sa 250 km above the base of the 
photosphere) . 

From the available data wc use only the SOLIS /VSM magnetogram for 24 October 
and the Hinode/SP magnetogram for 27 October to derive the relevant boundary 
conditions for force- free modeling, as outlined in Scction l^TTl We exclude the data for 28 
and 29 October because the active region was too close to the limb, where projections 
effects become prohibitive. We exclude the Hinode/SP data for 24-26 October because 
for these days most of the active region lies outside the Hinode/SP field of view. We 
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exclude the SOLIS/VSM data for 25 and 26 October because the magnetograms appear 
to contain significant systematic errors and are not suitable for force-free modeling. In 
particular the component of the magnetic field transverse to the line of sight shows a 
strongly preferred direction which we attribute to the result of systematic error. 

The lower panel of Figure [5] shows the times of the acceptable magnetogram ob- 
servations (solid vertical lines) compared with the flare history of the region. The 
magnetogram data requires ~ 10 minutes of integration time in the case of SOLIS/VSM 
and w 60 minutes of integration time in the case of Hinode/SP. The region first began 
flaring on 24 October producing several small GOES A-class and B-class, and the 
SOLIS/VSM magnetogram data are taken on this day. The flaring rate for the region 
peaked on 26 and 27 October and the Hinode/SP magnetogram data are obtained 
during this interval. 



3.1. SOLIS/VSM Data 

The National Solar Observatory SOLIS/VSM provides Milne-Eddington inverted vec- 
tor magnetograms, and 'Quicklook' data consisting of init ial estimates of the vector 
magn etic field produced prior to Milne-Eddington inversion ( Hennev. Keller, and Harvevl .[ 



20061 ). Both full disk magnetograms and cropped data centered on individual active 
regions are available onlin£0. The online data provide the magnetic field vector close 
to the photosphere i n terms of a line-of-sight compone nt and a component transverse 



to the line of sight (jHennev. Keller, and Harvevl . l2006f ). The transverse component is 



further decomposed into a magnitude and direction defined by the azimuthal angle. The 



i nvers ion is performed by the NSO and details are given in lHennev. Keller, and Harvev 



(120061) 



The magnetogram for 24 October is obtained at 16: 05 (UT) and is 236 x 236 pixels 
in size with a spatial resolution of 1.1" per pixel. The magnetogram covers an area 
of 267" X 267" corresponding to 195 x 196 Mm on the photosphere. The NSO also 
provides Quicklook data for the same day. Over regions of the photosphere where the 
magnetic field is weak there are gaps in the Milne-Eddington data where inversion in 
not performed due to a poor signal to noise ratio. We embed Quicklook data directly 
into the Milne-Eddington data to fill these gaps, following the advice provided by the 

NscEI. 

As explained in Section l^TTI force-free modeling is performed in a locally fiat Cartesian 
coordinate system defined on the photosphere, and we refer to this coordinate system 
as lielioplanar coordinates. The magnetogram is provided on a uniform grid on a plane 
perpendicular to the line of sight. We map the magnetogram to a helioplanar system 



using the transformations given in IVenkatakrishnan. Hagvard. and Hathawavl (jl988l) 



The helioplanar coordinate system is defined on a plane tangent to the solar surface at 
a point at the center of the active region. The axes of this system are oriented such that 
the X and y axes align with solar west and north respectively, and the z axis is normal 
to the photosphere. Owing to projection effects after mapping, the grid on which the 
magnetic field is defined is no longer regularly spaced. We regrid the magnetogram 
onto a larger (in number of grid points) uniform grid with the same spatial resolution 
(1.1" per pixel) as the original data. Since the active region is close to disk center on 
24 October the projection effects are small. 



See|http:/ /solis. nso.edu.' 

See the introduction to SOLIS/VSM data at |http://solis.nso.edu/vsm/vsm_imageJnfo3.html[ 
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The spectropolarimetrically derived azimuthal angle contains an inherent 180° am- 
biguity (jLandi Degl'Innocenti and Landolfil . |2004[ ) which must be resolved to determine 
the magnetic field components {B^, By and B^) needed to derive boundary conditions 
for force-free r nodeli ng. Several procedures exist for resolving the 180° ambiguity [see 



Metcalf et al. 



minimum energy 



(|2006l) and refer e nces t herein], and here we use the 

method developed bv iMetcalf (1994). This method determines the azimuthal angle 
by requiring the ambiguity-resolved field to simultaneously minimize |V • B| and |Jp 
over the magnetogram. The optimal magnetic field con figuration is found using the 



general optimization procedure of simulated annealing (jPress et all 119921) . We note 



that the NSO/SOL I S pro vides ambiguity-resolved data using a different resolution 
method ( Geor goulid . l2005[) . but we have chosen to perform this step again using the 
minimum energy procedure. 

The boundary conditions for force-free modeling are the component of the magnetic 
field normal to the photosphere, which corresponds to the B^ values in the regridded 
magnetogram, and the force-free parameter a on the photosphere derived from the 
values of B^ and By in the magnetogram using Equation ([5l. The partia l deriv atives 



in Equation are approximated using centered differencing (jPress et al\ . Il992l ) . The 
boundary values on a are set to zero at points where \Bz\ < 0.01 x maxdiJ^j) because 
the ratio Jz/B^ is inaccurate in regions where B^ is small. 

The final boundary conditions after projection and regridding are shown in Figures[3] 
andm Figure [3] shows the boundary conditions for B^ (upper panel), and the boundary 
conditions for (lower panel). Figure |4] shows the force- free parameter ao, and the 
signal to noise ratio \ao\/a for ao where a represents the approximate uncertainty in 
the ao value estimated using the B^ boundary value, as explained in Section 13.41 In 
both Figures [3] and 01 the data saturate at a particular threshold value which is below 
the maximum in the data, by which we mean all points above the threshold appear 
with the same color. This is done to provide a good contrast between large and small 
data values, and to diminish the appearance of outliers. For the panels showing and 
ao only points where the signal to noise ratio |ao|/cr is three standard deviations above 
the mean signal to noise ratio are represented, points where |ao|/cr = are excluded 
when computing the mean and standard deviation. This criteria is chosen because it 
is independent of the magnitude of cr which is arbitrary for the SOLIS/VSM data (see 
Section 13. 4p . Only the central portion of the data is shown in the two panels because 
the boundary data values (S^, Jz and ao) outside the central region are comparatively 
small and are not visible with the chosen scale. 

3.2. Hinode/SP Magnetogram Data 

The Hinode/SP vector magnetogram data for 27 October is derived from spectropolari- 
metric d ata from the Solar Optical Telescope/Spe ctroPolarimeter aboard the Hinode 



satellite (|Kosugi et all l2007l : iTsuneta et M . l2008l ). The 'fast -scan' map obtained at 



10:45(UT) is spatially sampled at 0.32", and covers a relatively small 360" x 164" 
field of view. The data are inverted using the High Altitude Observatory (HAO) Milne- 



ligi .. . 

Eddi ngton inversion code ( Skumanich and Litesl . 19871 : iLites and Skumanichl . ll990l : Lites et 
Il993l ). which has been modified for Hinode/SP data (revised code courtesy B. Lites, 
HAO/NCAR). Ambiguity r esolution is performed using an implementatior|§ of the 
minimum energy procedure ( Leka. Barnes, and Crouchl . l2009h . 

For the Hinode/SP data for 27 October uncertainties in the three components of the 
magnetic field {a^-, cry,az) are derived during the Milne-Eddington inversion. These are 



^Available at |http:/ /www. cora.nwra.com /AMBIG/I 
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Figure 3. The central part of part of the SOLIS/VSM vector magnetogram data for 24 October. 
The upper panel shows and the lower panel shows J^. The data are colored blue when positive, 
white when zero, and red when negative. Contours of are superimposed on the data, with contours 
separated by intervals of 300 G. The contours are black where B^ < and light gray where B^, > 0. 
The values of B^ (upper panel) have been saturated at ±944 G with a maximum value of 961 G. The 
values of Jz (lower panel) have been saturated at itlSmAm"'^ with a maximum value of 23mAm~'^. 
In the lower panel only points where \ao\/cr is three standard deviations above the mean are shown. 
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Figure 4. Boundary conditions on the force- free parameter obtained from the SOLIS/VSM vector 
magnetogram data for 24 October. The upper panel shows ao and the lower panel shows the signal to 
noise ratio |ao|/a". The data are colored blue when positive, white when zero, and red when negative. 
Contours of Bz are superimposed on the data, with the contours separated by 300 G. The contours 
are black where Bz < and light gray where Bz > 0. The values of ao (upper panel) have been 
saturated at ±lMm~^ with a maximum value of 3.7Mm~^. The values of jaol/f (lower panel) have 
been saturated at 50 with a maximum value of 67. In the upper panel only points where |q:o|/c is 
three standard deviations above the mean are shown. 
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based on the uncertainty in the fit of the synthetic Stokes profiles to the observed 
Stokes profiles. The uncertainties are then propagated through the projection, regrid- 
ding and ambiguity resolution to derive uncertainties for aQ. These uncertainties are 
a measure of how well the artificial Milne-Eddington profiles fit the observed profiles. 
They do not take into account the approximation of the Milne-Eddington atmosphere 
or systematic errors in the observations of the Stokes profiles. Therefore they may 
be considered lower bounds to the real uncertainties. Some parameter-dependence is 
intro duced into the me thod of simulated annealing by the choice of a 'cooling sched- 
ule' (jPress et al\ . Il992r ) . We note that some data points are particularly sensitive to 



the choice of schedule, but overall the large scale structure of the ambiguity-resolved 
magnetogram appears independent of the choice. Large uncertainties are assigned to 
points which depend strongly on the choice of cooling schedule to reflect this parameter 
dependence. 

The resulting vector magnetic components and uncertainties from the Hinode/SP 
data are re-mapped to a helioplanar 1" grid. The Hinode/SP field of view is small and 
only covers the main portion of the active region. As explained in Section 12.11 when 
performing force-free modeling it is important to isolate the region of interest from 
the boundaries to limit the influence of the boundaries on the solution. As with other 
studies involving Hinode/SP data (e.g. Wheatland and Lekall2011 : De Rosa et fflZ.|[2009l ) 
the field of view is expanded by adding line-of-sight data from the Miche lson Doppler 



Imager on the Solar and Heliospheric Observatory (jScherrer et all . 119951) around the 



Hinode/SP data. A map of Bz is constructed using a potential field matching the 
observed line-of-sight magnetic field values in the boundary, and this is also remapped 
onto a 1" grid. After coalignment between the two data sources, an apodizing function 
is applied at the boundaries of the t wo data sets to ensure a smo oth transition between 
data values. Th is procedure fol l ows Wheatland and Lekal (j201lr ). and improves upon a 
method used in lDe Rosa et ali (|2009l ). The final size of the boundary data is 440 x 300 
pixels, covering a region on the photosphere of size 320 x 220 Mm. 

The data are shown in Figures [5] and [6l The vertical magnetic field Bz (upper panel) 
and vertical current density Jz (lower panel) are shown in Figure [5] The presentation 
of the data follows the format used in Figures [3] and SI Approximate uncertainties a 
in the ao values are provided by the spectropolarimctric inversion, for the Hinode/SP 
data points. Figure [S] shows the force- free parameter ao (upper panel) and the signal- 
to- noise ratio |ao|/cr (lower panel). For the panels showing Jz and ao only points where 
the signal to noise ratio \ao\/a is three standard deviations above the mean signal to 
noise ratio are represented, points where laol/f = are excluded when computing the 
mean and standard deviation. Only the central portion of the data is shown in the two 
panels because the boundary data values {Bz, Jz and ao) outside the central region 
are comparatively small and are not visible with the chosen scale. 



3.3. Flux Balance 



The integral of Bz over the magnetogram defines the net magnetic fiux across the 
observed region of the photosphere. If the net fiux is zero then the magnetogram is said 
to be fiux balanced. Flux balance is not essential for the force-free modeling procedure 
used here, but if the flux is unbalanced a greater number of fleld lines leave the top and 
sides of the computational domain in order to satisfy V • B = globally. As discussed in 
Section 12.11 boundary conditions are not known on the top and side boundaries of the 
domain and the force-free code sets a = along fleld lines which leave the domain via 



these boundaries (jWheatlandl . 120071 ). The currents and therefore the free energy tend 
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Figure 5. The central part of the Hinode/SP vector magnctogram data for 27 October. The upper 
panel shows and the lower panel shows J^. The data are colored blue when positive, white when 
zero, and red when negative. Contours of are superimposed on the data, with contours separated 
by intervals of 1000 G. The contours are black where B^ < and light gray where Bz > 0. The values 
of Bz (upper panel) have been saturated at ±2448 G with a maximum value of 2388 G. The values of 
Jz (lower panel) have been saturated at itSOmAm"^ with a maximum value of 144mAm~^. In the 
lower panel only points where laol/o" is three standard deviations above the mean are shown. 



to be reduced by this procedure. The effect of the boundaries on the energy estimates 
is likely to be less pronounced if the magnctogram is flux balanced. 

A practical measure of flux balance is the relative flux balance, defined as 

f B,dA 

where the integrals are over the magnctogram area. The integrals in Equation ([7]) are 
evaluated numerically using the two-dimensional trapezoidal rule. The results are listed 
in the right-hand column in Table [TJ 

We consider a $b of less than 1% acceptable. The Hinode/SP magnctogram for 27 
October meets this requirement, but the SOLIS/VSM magnctogram for 24 October 
does not. To flux balance the SOLIS/VSM data a buffer region of constant is 
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Figure 6. Boundary conditions on the forcc-frcc parameter obtained from the Hinode/SP vector 
magnetogram data for 27 October. The upper panel shows ao and the lower panel shows the signal to 
noise ratio \ao\/cr. The data are colored blue when positive, white when zero, and red when negative. 
Contours of Bz are superimposed on the data, with the contours separated by 1000 G. The contours 
are black where Bz < and light gray where Bz > 0. The values of «o (upper panel) have been 
saturated at ±0.15 Mm"-*^ with a maximum value of 1.9 Mm~^. The values of \ao\/cr (lower panel) 
have been saturated at five with a maximum value of 5.4. In the upper panel only points where \ao\/a 
is three standard deviations above the mean are shown. 



added around the border of the magnetogram. The magnitude of the field in the buffer 
region satisfies < O.Of x maxdS^j), and hence the buffer region is considered weak 
field by our previous definition (see Section 13. fp . In the buffer the magnitude of 
is 9.3 gauss (G) (the weak field threshold is 9.6 G), and the value of uq is set to zero 
everywhere in the buffer. The total size of the magnetogram with the buffer added is 
270 X 270 grid points, and the resolution of the data remains the same (l.f " per pixel). 
The final size of the region including the flux buffer is 222 x 222 Mm. 

3.4. Phenomenological Uncertainties for SOLIS/VSM Data 

Uncertainties are provided with the Hinode/SP data based on spectropolarimetric 
inversion of the data (as explained in Section l3.2p . but are not available for the SO- 
LIS/VSM data for 24 October. Instead we make a simple estimate for the uncertainties 
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based on the magnitude of at each boundary point. Specificany we introduce a 
phenomenological uncertainty defined by 



{0.01 for the flux buffer region; 
B,,,/\B,{x,y)\ where \B,/B^\ > 0.01 (8) 
100 where \BJBJ\ < 0.01 

Equation ([5]) is written in non-dimensional units. The magnetic field is scaled by _Bi„ = 
niaxdS^I), the maximum value of \Bz\, and a is scaled by L^^ the transverse length of 
the magnetogram. Points in the flux buffer are assigned a small constant uncertainty, 
which ensures ao = for this region during the self-consistency procedure. A relatively 
large uncertainty {Lxcr ~ 100, in non-dimensional units) is assigned to points in weak 
field regions where ao has been set to zero (see Section [3TT|) . 

The phenomenological scaling a{x, y) ~ l/\Bz{x, y)\ for the SOLIS/VSM data region 
is based on inspection of the inversion-based uncertainties for the Hinode/SP data for 
27 October. 

Figure [7] shows P{a, where P{a, \Bz\) is the distribution of the inversion- 

based uncertainties for the Hinode/SP data for 27 October over field strength and 
uncertainty in ao i.e. the fraction of points with field value \Bz\ and uncertainty a. 
We have applied a nonlinear scaling to the distribution so that strong field regions 
with small values of P{a, \Bz\) are visible. The solid line shows a = 100 Mm"-'^/!^^! 
with Bz in gauss, to illustrate cr l/l^zl dependence for comparison. The choice of 
the factor 100 Mm~^ is arbitrary. The figure shows that the simple phenomenological 
scaling implied by Equation is accurate for a majority of data points, although their 
is notable departure for large values of \Bz\. 

It should be noted that, according to Equation the exact magnitude of the un- 
certainties is not important for the self-consistency averaging procedure [see discussion 
following Equation 1^]. Only the relative size of uncertainties at different points is 
important. 

Although, due to the departure from the a ^ l/|Sz| scaling for large values of \Bz\ 
we expect that the uncertainties in strong field regions are underestimated relative to 
weak field regions. 



3.5. Force-free Modeling Using SOLIS/VSM and Hmode/SP Data 



Force-free solutions are calculated using the code of IWheatlandl ()2007l) and the self 



consistency procedure of Iwheatland and Regnierl ( 2009t ). The boundary conditions 



derived from the SOLIS/VSM and Hinode/SP vector magnetogram data described 
in Sections 13.11 and 13.21 are used. The initial magnetic field for the Grad-Rubin it- 
eration is a potential field calculated using a Fourier solution to Poisson's equation. 
The computational domain used is a uniform three dimensional Cartesian grid with 
transverse dimensions (i.e. the x and y dimensions) matching the dimensions of the 
magnetogram which forms the lower boundary at z = 0. The vertical dimension (i.e. 
z) is 128 points high for the SOLIS/VSM data and 200 points high for the Hinode/SP 
data. The resulting dimensions of the grids are 270 x 270 x 128 points for 24 October, 
and 440 x 300 x 200 points for 27 October. The grid spacing is uniform in each direction. 
We summarize the size, maximum magnetic field, and flux balance of the data in Table 

m 
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Figure 7. The relationship between the uncertainties on the boundary values of the force-free param- 
eter derived from the spectropolarimetric inversion and \B~\ for the Hinode/SP data for 27 October. 
The solid line is to illustrate cr ~ l/IB^I dependence and the position of the line has been chosen 
arbitrarily. There is a noticeable departure from the cr ~ 1 / 1 B, \ scaling for large values of \Bz\. 



Table 1. Summary of magnctogram data/computational data domain sizes, field strength parameters,! 
and relative flux balance (sec Equation lO) for the computational models of AR 11029. | 





Data source 


Date 


X Ny X TVz 


Lx X Ly X Lz 

(10* m) 


Bm = max( Bz \ ) 
(G) 


(%) 


SOLIS/VSM 
SOLIS/VSM (buffered) 
Hinode/SP 


Oct. 24 
Oct. 24 
Oct. 27 


236 X 236 X 128 
270 X 270 X 128 
440 X 300 X 200 


2.0 X 2.0 X 1.1 
2.2 X 2.2 X 1.1 
3.2 X 2.2 X 1.5 


962 
962 
2448 


11 

0.02 

0.4 



4. Results 

This section presents the results of the force-free modehng of active region AR 11029. 
Section Wl] discusses the results when the force-free code is applied to the SOLIS/VSM 
data for 24 October to construct P and TV solutions without using the self-consistency 
procedure. Section discusses the application of the self-consistency procedure to the 
same data. Section 14.31 discusses the application of the self-consistency procedure to 
the Hinode/SP data for 27 October. 

4.1. Construction of P and N Solutions Directly from Magnetogram Data for 24 
October 

We construct P and solutions for AR 11029 directly from the magnctogram data for 
24 October. In this section wc do not apply the self-consistency procedure. To construct 
the P and N solutions 50 Grad-Rubin iterations are used beginning with a potential 
field constructed from B^. 
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Figure 8. Results for the construction of P and N solutions directly from the SOLIS/VSM magne- 
togram data for 24 October. The total magnetic energy is shown (as a fraction of the energy of the 
potential field) for the P solution (plus signs) and the N solution (diamonds) over 50 Grad-Rubin 
iterations starting from a potential field. The oscillations in the energy occur because of large localized 
currents in the boundary data. 



In application to the magnetogram data the Grad-Rubin procedure does not strictly 
converge for the P and N solutions. The energy of the magnetic field at each iteration 
provides a way to examine the convergence. Figure [5] shows the magnetic energy, E, in 
units of the energy of the potential field, Eq, over successive Grad-Rubin iterations for 
both the P and N solutions. The energy of the two solutions behaves similarly. Initially 
the total energy increases beyond the energy of the potential field and then begins 
oscillating after about 15 iterations for the P solution and after about ten iterations 
for the N solution. This oscillatory behavior is typical in situations where the Grad- 
Rubin method does not strictly converge, and may be attributed to large currents at 
certain points in the boundary. This is confirmed by inspecting visualizations of field 
lines at successive Grad-Rubin iterations. For the P and N solutions for the 24 October 
data we observe that the longest field lines which stretch high into the computational 
volume (and carry significant currents) change over successive iterations. The field in 
the lower part of the volume is static, and the field in the middle of the volume displays 
small oscillations. 

We estimate the energy of the P and N solutions by averaging the energy over several 
Grad-Rubin iterations. After 15 Grad-Rubin iterations the P solutions stops changing 
systematically and begins oscillating so we estimate the energy of the P solution by 
taking the average energy of the field over successive Grad-Rubin iterations starting 
from iteration 15 [i.e. wc take Ep = mean(i?i5, i^ig, i?5o), where Ei is the energy 
of the P solution after i Grad-Rubin iterations]. We do the same for the solution, 
except we start the average from iteration ten [i.e. En = mean(i?io, £'50)]- We also 
compute the free energy E[ = E — Eq, for E = Ep and E = E^. The results are listed 
in Table H 

The upper panel of Figure |9] shows a visualization of the magnetic field after 20 
Grad-Rubin iterations (the lower panel is discussed in Section IT^ . The magnetic field 
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Table 2. Energy estimates for the P and A'^ so- 
lutions constructed directly from the SOLIS/VSM 
magnetogram data for AR 11029 on 24 October. 



Solution E Eo Ef = E - Eg 

(10^1 erg) (lO^'ierg) (10^" erg) 

P 9.98 9.62 3.6 

N 9.77 9.62 1.5 



lines for the P solution (blue curves) and the N solution (red curves) are displayed 
in the upper panel. The fields are viewed looking down on the photosphere, which is 
colored blue where > and red where < 0. The panel shows some qualitative 
difference between the P and N solutions. 

The visualization of the field lines also reveals the extent to which the missing 
information on the boundaries influences the solution. It can be seen from the upper 
panel of Figure that both solutions have field lines which cross the side boundaries. 
Many field lines connected to the trailing positive spot leave the box. However, a 
significant part of the active region field involves closed field lines. It is unlikely the 
boundary has a large infiuence on the final energy. 

The results in this section illustrate difficulties associated with obtaining estimates 
for the total and free energy using the nonlinear force-free model. The Grad-Rubin 
iteration does not strictly converge, due to strong currents, which introduces a degree 
of arbitrariness in the energy estimate since the energy depends on the choice of when 
to stop the iteration process. Compared with the total energy the oscillations in the 
total energy are small, but these oscillations in the total energy translate into large 
oscillations in the free energy, because E ^ E{. The force- free energy estimates may 
be regarded as approximate (order of magnitude) estimates of the free energy. Based 
on the results in Table [2] we estimate that the free energy is « 3 x 10'^° erg. 

4.2. Construction of Self-consistent Solutions for 24 October 

The self-consistency procedure is applied to the SOLIS/VSM vector magnetogram 
boundary data for 24 October. Initially, we use ten self-consistency cycles with 20 
Grad-Rubin iterations per cycle (P and N solutions are constructed using 20 Grad- 
Rubin iterations, and this is repeated for ten self-consistency cycles, with the ao-valucs 
updated using Equation ([5]) at the end of each cycle). 

The energy of the P and N solution at each self-consistency cycle may be used to 
monitor the convergence of the procedure. At the end of each cycle (after 20 Grad- 
Rubin iterations for both P and N solutions), the total magnetic energy, E, for the 
solutions is calculated. Figure [TU] shows 

the evolution oi E/Eq with cycle for the P and N solutions. Initially the two energies 
are different (as discussed in Section HTTI and shown in Tabled but after three cycles 
the energy approaches a single value. After ten cycles E/Eq = 1.0047 for the P solution 
and E/Eq = 1.0044 for the solution which is lower than the energy of either the 
initial P or N solutions. The small difference between the energies of the P and A^ 
solutions translates into a difference of « 7% in the free energy. In dimensional units 
the free energy is about 4 x 10^^ erg. The total energy of the final self-consistent solution 
E, the energy of the potential field Eq, and the free energy of the final self-consistent 
solution E{ = E — Eq are listed in Table 12] (the first two rows of this table are relevant; 
the second two rows are discussed below). 
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Figure 9. Nonlinear force-free models of the coronal magnetic field constructed from the SOLIS/VSM 
vector magnctogram data for 24 October. The P solution (blue field lines) and N solution (red field 
lines) arc shown viewed looking down on the photosphere. The upper panel shows the two solutions 
constructed from the original vector magnetogram boundary data (with a flux buffer added) using 20 
Grad-Rubin iterations. The lower panel shows the solutions after ten self-consistency cycles with 20 
Grad-Rubin iterations per cycle. The photosphere is colored blue where Bz > and red where Bz < 0. 
This figure illustrates the achievement of self-consistency: the P and A'^ solutions in the upper panel 
are dissimilar but in the lower panel are very similar. 



The lower panel of Figure shows a visualization of the magnetic field of the self- 
consistent solutions. The magnetic field lines of the P and TV solutions after ten self- 
consistency cycles are shown by the blue and red curves respectively. In contrast to 
the top panel of Figure IHl the two sets of field lines are very similar. Small differences 
are seen for field lines extending high into the computational volume. The agreement 
between the two sets of field lines confirms that a single self-consistent solution is 
obtained. 
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Figure 10. Results for the application of the self-consistency procedure to the SOLIS/VSM vector 
magnetogram data for 24 October. The energies of the P and solutions (as a fraction of the energy 
of the potential field) are shown over the ten self-consistency cycles. For each cycle 20 Grad-Rubin 
iterations are used to construct the P and N solutions. The energy of the P solution is indicated by 
plus signs and a dashed line, and that of the TV solution by diamonds and a solid line. 

Table 3. Energy estimates for the self-consistent solutions constructed from the 
SOLIS/VSM magnetogram data for AR 11029 on 24 October. 



Grad-Rubin iterations / cycle 


Solution 


E 


Eo 


Ei = 






(10^1 erg ) 


(10^1 erg) 


(1029 


20 


P 


9.669 


9.62 


4.52 




N 


9.667 


9.62 


4.23 


35 


P 


9.666 


9.62 


4.20 




N 


9.667 


9.62 


4.23 



The oscillations of the energies of the initial P and N solutions (see Section 14. ip 
means that the magnetic field used to re-map ao at the end of the early self-consistency 
cycle depends somewhat on the number of Grad-Rubin iterations used to construct the 
P and N solutions. If a different number of Grad-Rubin iterations is chosen the P 
and N solutions have different field lines, because of the oscillations seen in Figure [51 
Hence the mapping of field lines is different and the application of Equation ^ gives a 
different result. Correspondingly the final energy of the self-consistency procedure de- 
pends somewhat on the number of Grad-Rubin iterations used. To test this dependence 
we repeat the calculation with 35 Grad-Rubin iterations for the P and N solutions, 
instead of 20. We choose this value because it corresponds to an iteration at which 
the difference in energy between the P and TV solutions constructed directly from the 
vector magnetogram data is large (sec Figure [5]). Once again ten self-consistency cycles 
are applied. The final energy and free energy in this case arc also given in Tabic [3] 
in the third and fourth rows. The variation in the energies for the two self-consistent 
solutions with different numbers of Grad-Rubin iterations per self-consistency cycle 
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is sufficiently small to justify an order of magnitude (or better) estimate for the free 
energy of the region based on the results. 

The self-consistency procedure modifies ao according to the procedure explained in 
Section 12.21 Figure [TT] shows the vertical current density after ten self-consistency 
cycles for the P solution (upper panel) and for the N solution (lower panel). The self- 
consistent values of Jz for the P and N solutions are very similar, confirming that a 
self-consistent force-free solution is found. This figure should be compared with the 
lower panel of Figure [3j which shows for the original magnetogram data. Both 
Figures [51 and [TTI cover the same region of the photosphere. 

The currents in the lower boundary are significantly reduced as a result of the 
self-consistency procedure (the maximum currents are « 5 mAm^^ compared with 
sa 15mAm^^ in the original data). This is in part due to the averaging implied by 
Equation ((5]). However, it also represents the preservation of smaller but more certain 
values and the loss of larger but less certain values (based on the phenomenological 
uncertainty rule a ^ \Bz\~^). Specifically the values of ao are typically smaller in 
strong field regions, which have a smaller uncertainty, so these smaller values are better 
preserved in the application of Equation ([6]) (compared to large values of ao in weak field 
regions). There are substantial areas with ag = in the final self-consistent boundary 
data. These patches correspond to open field regions where a is set to zero by the 
handling of the side and top boundary conditions, as explained in Section \2A\ 

4.3. Construction of Self-consistent Solutions for 27 October 

The self-consistency procedure is applied to the Hinode/SP magnetogram data for 27 
October. The initial P and N solutions again oscillate in energy. However, for the 
Hinode/SP data the oscillations are so severe that a meaningful estimate of the energy 
and free energy cannot be made. The large oscillations are due to large boundary 
values of ao. Many of these large values have a large associated uncertainty, and hence 
may be spurious, but the results of the calculation are infiuenced by these values. It is 
necessary to apply the self-consistency procedure which accounts for uncertainties in 
ap and tends to remove values with large uncertainties in a systematic fashion. 

The self-consistency procedure is applied to the Hinode/SP data for 27 October 
using two separate runs with a different number of Grad-Rubin iterations per cycle 
for each run, similar to the procedure in Section 14.21 with the SOLIS /VSM data. For 
the first run we use 20 Grad-Rubin iterations per self-consistency cycle, and for the 
second run we use 30 Grad-Rubin iterations per self-consistency cycle. We apply ten 
self-consistency cycles in each case. For both runs the final energies of the P and TV 
solution converge to a single value after ten self-consistency cycles. Figure [T^] shows 
the energy of the P solution (plus signs) and N solution (diamonds) in units of the 
potential field energy at the end of each self-consistency cycle over the ten cycles with 
20 Grad-Rubin iterations per cycle. The results using 30 Grad-Rubin iterations are 
similar to those using 20 iterations per cycle. In both cases the energies of the P and N 
solution initially differ greatly, but after ten self-consistency cycles the energy difference 
is reduced to « 1%. 

The total energy, free energy of the P and N solutions, and the energy of the 
potential field for the self-consistency calculation arc summarized in Table 01 

The total and free energy estimates differ for the two runs with different numbers 
of Grad-Rubin iterations, but the difference is sufficiently small to allow an order of 
magnitude estimate for the energy. The total energy is large (of order lO'^'^ erg, which is 
ten times larger than was found for the region on 24 October). The variation in the free 
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Figure 11. Results for the application of the self-consistency procedure to the SOLIS/VSM vector 
magnotogram data for 24 October. The vertical electric current density Jz over part of the boundary is 
shown for the P solution (upper panel), and for the N solution (lower panel) after ten self-consistency 
cycles. The values of Jz are saturated at ±5mAm~^, and only points where \ao/a\ is three standard 
deviations above the mean are shown. The contours are the same as Figure [3] The distributions of 
Jz for the P and A'^ solutions after the self-consistency procedure are very similar, indicating a single 
self-consistent solution is found. 
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Figure 12. Results for the application of the self-consistency procedure to the Hinode/SP vector 
magnetogram data for 27 October. The energies of the P and N solutions (as a fraction of the energy 
of the potential field) are shown over the ten self-consistency cycles. For each cycle 20 Grad-Rubin 
iterations are used to construct the P and N solutions. The energy of the P solution is indicated by 
plus signs and a dashed line, and that of the TV solution by diamonds and a solid line. 

Table 4. Energy estimates for the self-consistent solutions constructed from the 
Hinode/SP magnetogram data for AR 11029 on 27 October. 



Grad-Rubin iterations / cycle 


Solution 


E 


Eo 


Ei = 


E-Eo 






(10^3 erg ) 


(10^3 erg) 


(1031 


erg) 


20 


P 


1.769 


1.707 


6.16 






N 


1.772 


1.707 


6.50 




30 


P 


1.787 


1.707 


7.94 






N 


1.791 


1.707 


8.35 





energy is correspondingly large, but we can still place an upper bound of « 10^^ erg on 
the free energy of the active region based on the results in Table HI 

Figure [T3] shows a visualization of the magnetic field of the self-consistent P solution 
(blue curves) and N solution (red curves) for the run using 20 Grad-Rubin iterations 
per cycle. The two sets of field lines are similar, although some small differences are 
visible. 

The self-consistency procedure modifies ag according to the procedure explained in 
Section [321 Figure [HI shows the vertical current density Jz after ten self-consistency 
cycles for the P solution (upper panel) and for the N solution (lower panel). The 
results are shown for the run using 20 Grad-Rubin iterations per cycle, but the results 
using 30 Grad-Rubin iterations per cycle are similar. The self-consistent values of Jz 
for the P and N solutions are very similar, confirming that a self-consistent solution 
is found. The self-consistent values of Jz are significantly reduced by comparison with 
the values of Jz in the original magnetogram (lower panel of Figure [S|) . This is due 
to the averaging implicit in the procedure, and the removal by the procedure of large 
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Figure 13. Nonlinear force-free models of the coronal magnetic field constructed from the Hinode/SP 
vector magnetogram data for 27 October. The self-consistent P solution (blue field lines) and N solu- 
tion (red field lines) are shown for the calculation using ten self-consistency cycles with 20 Grad-Rubin 
iterations per cycle. The region is view looking down on the photosphere which is colored blue where 
Bz > and red where < 0. 



values with large associated uncertainties (see discussion in Section [^T^ . The final self- 
consistent boundary data contains regions with ao = corresponding to open field 
regions where a has been set to zero by the handling of the side and top boundary 
conditions (sec Section ■ 



5. Summary and Discussion of Results 

Estimates of the total energy and free energy of the coronal magnetic field of active 
region AR 11029 are calculated using two vector magnetograms. The first magne- 
togram is based on data for 24 October from the SOLIS Vcctor-SpectroMagnetograph 
(SOLIS/VSM), and the second is based on data for 27 October obtained from the 
Hinode Solar Optical Telescope SpectroPolarimeter (Hinode/SP). The estimates use 
nonlinear force- free solutions constructed from the vector magnetogram boundary data. 
The magnetogram data provide two solutions to the nonlinear force-free model (the 
P and N solutions, denoting the choice of polarity in the boundary conditions used). 
The model does not have a unique solution because the boundary data are inconsistent 
with the model. In each case 'self-consistent' solutions ( Wheatland and Regnier . 20091 ) 



arc constructed taking into account uncertainties in the boundary data. For the 24 Oc- 
tober data the uncertainties are estimated phenomenologically, and for the 27 October 
data the uncertainties are obtained from the inversion process used to construct the 
magnetogram. 
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Figure 14. Results for the application of the self-consistency procedure to the Hinode/SP vector 
magnetogram data for 27 October. The vertical electric current density Jz over part of the boundary 
is shown for the P solution (upper panel), and for the A'^ solution (lower panel) after ten self-consistency 
cycles. The values of Jz are saturated at ±30mAm~'^, and only points where |ao/'''l is three standard 
deviations above the mean are shown. The contours are the same as Figure [S] The distributions of 
Jz for the P and N solutions after the self-consistency procedure are similar, indicating a single 
self- consistent solution is found. 



The initial P and N solutions calculated from the vector magnetogram data for 24 
October have different energies and different field line structures. These solutions are 
themselves not uniquely defined because the Grad-Rubin iteration procedure used to 
construct the solutions does not strictly converge. The energy of the solutions oscillate 
with successive Grad-Rubin iterations owing to large currents in the boundary data. 
The amplitude of the oscillations in the total energy is relatively small, but the free 
energy is a small fraction of the total energy so this represents large amplitude oscilla- 
tions in the free energy. The oscillations limit the accuracy of the estimation of the free 
energy, but an order of magnitude estimate is possible. The free energy of AR 11029 
on 24 October is estimated from the initial P and N solutions to be w 3 x lO'^" erg. 

The self-consistency procedure is applied to the 24 October data. The energy of the 
P and N solutions agrees closely after ten self-consistency cycles, and visual inspection 
of the field lines of the P and N solutions confirms that a self-consistent solution is 
obtained. The final energy of the self-consistent solution depends only weakly on the 
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number of Grad-Rubin iterations per self-consistency cycle allowing a reliable estimate 
of the free energy. The free energy based on the self-consistent force-free solutions for 
24 October is « 4 x 10^^ erg. This value is roughly an order of magnitude smaller than 
the value obtained from the initial P and N solutions. 

For the 27 October Hinode/SP data, the initial P and N solutions do not permit 
a meaningful estimate of the free energy, because of large amplitude oscillations in 
the energy during the Grad-Rubin procedure. The misbehavior is attributed to large 
currents in the boundary data, and many of these values also have large associated 
uncertainties. The large values prevent convergence of the Grad-Rubin iteration, and 
hence the result is completely determined by the spurious boundary data. We believe 
the self-consistency procedure which accounts for these large uncertainties provides a 
more reliable result. The self-consistency procedure is applied to the Hinode/SP data 
leading to the estimate that the free energy of the region is w 7 x 10'^^ erg on 27 October. 

The results present ed here illustrate the practical difficulties with force-free modeling 
( De Rosa et al. . 20091) . The models constructed directly from the vector magnctogram 



data do not provide a unique energy estimate so we use the self-consistency procedure. 
The accuracy of the energy estimates obtained using the self-consistency procedure has 
not been tested in any independent way, and it is possible to identify aspects of the 
modeling which may lead to underestimation of the energy. A number of factors tend 
to reduce the currents in the computational volume and in turn tend to reduce the free 
energy. Missing boundary data on the top and sides of the computational volume is 
handled by setting a = along field lines which leave the top and sides of the volume. 
This step removes currents from the volume. The self-consistency procedure also tends 
to reduce the ao values on the lower boundary. This is due to three reasons. First, 
large values of ag which have a large uncertainty are averaged out by the weighted 
average in the self-consistency procedure. Second, regions with ao ^ connected by 
field lines to regions in the flux buffer with ao ~ tend to get averaged to ao = 
over successive self-consistency cycles because points in the flux buffer are assigned 
small uncertainties. Third, the boundary data have many small regions of ao with 
opposite sign and averaging between these regions will lead to a net decrease in ao over 
successive self-consistency cycles. 

There is also considerable uncertainty associated with the boundary magnetogram 
data. The fleld values are based on inversion of spectropolarimetric measurements with 
flnite spatial, depth, and temporal resolution, and the inversions in this case assume a 
Milnc-Eddington atmosph ere and the analytic form of the Stokes sp ectra described by 
Unno and Rachkofsky (e.o. lLandi Degl'Innocenti and Landolfil[200l . The active region 



considered here, NOAA AR 11029, evolved rapidly and was highly dynamic, producing 
Stokes spectra which were strongly Doppler shifted and multi-lobed over much of 
the region. The latter indicates that the observed spectra arc not consistent with a 
Milne-Eddington atmosphere, hence the returned field values may not be appropriate 
estimates of the underlying field. The inversion based uncertainties (explained in Sec- 
tion [32]) for the Hinode/SP data reflect this somewhat but arc likely underestimates. 
The inversions for the SOLIS/VSM data can be assumed to harbor the same concerns. 
There is also uncertainty associated w ith the resolution of the 180-degree ambiguity 
(jMetcaff et ad bood lleka et alihOO^ . 

Given these caveats, it is nevertheless interesting to compare the energy of the 
largest flares produced by AR 11029 with the free energy estimates derived for the 
region, although only a fraction of the free energy is expected to be released in any 
single flare. To make this comparison we flrst establish an approximate correlation 
between GOES soft X-ray peak flux and flare energy. One measure of flare energy 
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Figure 15. The relationship between RHESSI energy estimates for the non-thermal energy of a 
sample of flares from the literature (see Section^ and the GOES peak flux of the flares. The plot also 
indicates the GOES classes. The vertical dashed line shows the peak flux of a GOES C2.2 class flare 
(peak flux of 2.2 X 10^'' Wm-^)^ the largest soft X-ray flare produced by AR 11029. 



is the ener gy of non-thermal elec trons, which may be inferred from hard X-ray ob- 



servations ()Lin and Hudsonl . 119711) . We consider a sample from the literature of flare 
energies estimated from hard X-ray observations, restricting attention (for uniformity) 
to observations made by the Re uven Ramaty High-Energy Solar Spectroscopic Im- 



ager (RHESSI) ijLin et all 120021 ) . These flare energies may be compared with the 
observed GOES peak flux of each event to establis h an energy-peak flux s caling . 
Figure [TSl sh ows the flare en e rgy estimates made by Sain t-Hilairc and Bcnz) (|200^, 
Saiiit-Hilairc and Bend (|2005[ ). [Maltagliati. Falchi. and Te"ria ca (2006), Holman (20051) 
and Emslie et all ( 2004 ) for 14 different RHESSI hard X-ray flares, plotted against the 
observed GOES soft X-ray peak flux of each flare. The plot shows that, to within an 
order of magnitude, a GOES B-class event has an energy of 10^® erg, a C-class event 
has energy an energy of 10^^ erg, an M-class event has an energy of 10'^° erg, and an 
X-class event has energy above 10"^^ erg. This provides a useful rule-of-thumb scaling 
law for flare energy estimates from GOES classification. These energy estimates only 
take into account the energy of non-thermal electrons producing hard X-rays in the 
flare. If an event is eruptive then a significan t amount of energy can be carried off by 
the mass ejection ( Gopalswamv et all . l2009l ). 
many small eruptions during its lifetim43- 

On 24 October AR 11029 produced four B class flares, the largest of which is a B4.9 
class event with a peak soft X-ray flux of 4.9 x 10^^ Wm^^. The energy of each of these 
flares is estimated to be of the order 10^^ erg, which is consistent with the upper bound 
of 4 X 10^^ erg derived from the self-consistent force-free modeling and the upper bound 
of 3 X lO'^^erg obtained from the initial P and N solutions (Tables [2] and [3]). 



and active region AR 11029 produced 



''See the LASCO CME catalog [http://cdaw.gsfc. nasa.gov/CME_list7] 
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On 27 October AR 11029 produced 23 soft X-ray flares. The largest flare produced 
by AR 11029 was a GOES class C2.2 which occurred at 00: 24 UT on 28 October. The 
estimate from the modeling for the free energy is w 7 x lO''^ erg (Tabled]), which is 
comparable to an X-class flare. 

Our estimates of the free energy of active region AR 11029 indicate that the region 
had sufficient energy to produce an M-class or X-class flare on 27 October, but not on 24 
October. It is interesting that the estimate for the energy on 27 October is substantially 
larger than that for 24 October. The photosphcric extent of the region grew rapidl y 
during the interval 24 October to 27 October as seen in Figure [Tland lWheatlandl (l2010l) . 
and it appears that the coronal magnetic fleld increased in size and energy accordingly. 
The magnetic flux of the region increased significantly between 24 and 27 October. 
On 24 October the magnetic flux was « 5 x 10^^ Mx over the positive polarity and 
« 4 X 10^^ Mx over the negative polarity. On 27 October the flux is « 30 x 10^^ Mx 
over both polarities, an increase of a factor of 6-7. 



6. Conclusion 



Solar active region AR 11029 was a highly flare productive sunspot region which 
emerged on the Sun in late October 2009. The region produced 73 soft X-ray flares 
as recorded in the event list for the Geostationary Operational Environmental Satel- 
lites [GOES) over a period of about a week (24 October to 1 November). Statistical 
analysis of these flares sho wed evidence for departure from the standard power-law 
frequency-size distribution ( Wheatlandl . 2010l ). Specifically, an absence of large events 
was indicated. This was conjectured to be d ue to the s mall active region having insuf- 
ficient magnetic energy to power large flares ([Wheatland . 2010). The aim of this paper 
is to test this hypothesis by estimating the free energy (i.e. the energy available for 
flaring) using nonlinear force-free modeling of the coronal magnetic field of the active 
region from vector magnetogram data. 

Estimates of the magnetic free energy of the coronal magnetic field of AR 11029 are 
presented using vector magnetogram data taken by the Vector-SpectroMagnetograph 
at the Solar Observatory's Synoptic Long term Investigations of the Sun facility (SO- 
LIS/VSM), and the Solar Optical Telescope SpectroPolarimeter on the Hinode satellite 
{Hinode/SP). Our data set consists of a SOLIS/VSM magnetogram taken on 24 October 
2009 and a Hinode/SP magnetogram taken on 27 October 2009. Force-free so l utions 
are constructed from the data using the Grad-Rubin code due to IWheatlandl (l2007h 
and the self-consistency procedure due to Wheatland and Regnier ( 2009t ). The proce- 
dure address the p roblem of the inconsistency between the force-free model and the 
photosphcric data ( Metcalf et al. 1995t De Rosa et al. 2009fl . 

For the two vector magnetograms free energy estimates are obtained. For 24 October 
we estimate a free energy of « 4 x 10^^ erg, which is roughly the size of a GOES C-class 
or small GOES M-class flare. For 27 October we estimate a free energy « 7 x 10'^^ erg, 
which is roughly the size of a GOES X-class flare. Although free-energy estimates made 
using non-linear force-free modeling are subject to signiflcant uncertainty we believe 
that the order of magnitude of these estimates is correct. 

The results in this paper do not support the hypothesis that there is an observable 
upper limit on the sizes of flares produced by small active regions due to insufficient 
free energy in the active region magnetic field. The modeling re-opens the interesting 
question of why this region did not produce any large flares given that it had sufficient 
energy to do so. 
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Appendix 

Glossary of terms. 

• Model: the nonlinear force-free model. 

• Solution: the magnetic field B = B(r) and the force-free parameter a = a{r) 
(where r is the position vector) which solve Equations ([3]) and in a domain 
with boundary conditions for and a (over a single polarity of B^) prescribed 
on the boundary of the domain. 

• P Solution: a solution where boundary conditions for a are prescribed at boundary 
points where B^ > 0. 

• N Solution: a solution where boundary conditions for a arc prescribed at boundary 
points where B^ < 0- 

• Iteration: one step in a procedure, e.g. one Grad-Rubin iteration is going from 
k to fc -|- 1 i n the Grad-Rubin or current-field iteration procedure, as defined in 
Wheatlandl |20o3). 



Self-consistency cycle: construction of P and Grad-Rubin solutions (involving 
a certain number of Gra d-Rubin iterations in each cas e) from given B^ and 
boundary conditions (see Wheatland and Regnier 2009I ). 



Fixed point /convergence: a fixed point in the Grad-Rubin proce dure is defined by 
BW = Bl'^+il for iterations k and k + 1 (see IWheatlmdIboOTt ). This represents 



convergence of the procedure (when met for all points in a computational volume) . 
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